A heuristic approach to the weakly interacting Bose gas 
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Some thermodynamic properties of weakly interacting Bose systems are derived from dimensional 
and heuristic arguments and thermodynamic relations, without resorting to statistical mechanics. 
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I. INTRODUCTION. 

In a previous article 1 I have derived the thermody- 
namic properties of ideal quantum gases solely from di- 
mensional arguments, the Pauli principle and thermody- 
namic relations, without resorting to statistical mechan- 
ics. In the present publication I will extend this approach 
to weakly interacting bosons. 

Starting point is the Hamiltonian of interacting bosons 
in d dimensions 



n = J d d x¥( 



2m 



(1) 



Here \I/ T (x) and ^(x) are the creation and annihilation 
operators of the Bose field in d-dimensions. m denotes 
the mass of the Bose particles, vo > is the strength 
of the interaction potential V(x) = t>o<5(x) of the bosons. 
Main reason for choosing here a strictly local interaction 
is to keep the number of parameters small, which only 
makes a dimensional analysis possible. However, for low 
energy phenomena this approximation is not unreason- 
able. 

It is convenient to express vq in terms of a scattering 
length b which follows from dimensional analysis: 



mvo 



l/(d-2) 



d>2. 



(2) 



Sd is a prefactor with S3 = 4ir. In d < 2, where ([2]) 
diverges with vo — > 0, the scattering length behaves as 
b ~ q^ 1 where q denotes the transferred momentum since 
there is always a bound stated. However, throughout this 
paper we will stay in d > 2 dimensions. 

For vo = the state of the system is completely de- 
scribed by the occupation numbers n p of the single par- 
ticle states which are characterized by their momentum 
p and energy eo(p) = p 2 /(2m). Since the particles of 
the ideal Bose gas do not occupy any volume and since 
there is no Pauli principle for bosons, in the ground state 
all particles can have the same position. One conclusion 
one can draw from this is that the pressure cannot de- 
pend on the volume^. There are no collective excitations. 
This situation is completely changed in the presence of 



even a very weak repulsive interaction. The particles 
have now a finite cross-section and occupy a finite vol- 
ume. Local density fluctuations will propagate because 
of the repulsive forces between bosons. This behavior re- 
sults in the occurrence of a branch of sound waves i.e. 
acoustic phonons. The existence of acoustic phonons in 
the weakly interacting system will be the main assump- 
tion of the present considerations. 

In the hydrodynamic limit, i.e. for frequencies ui <C 
T reiax wnere T reiax is b, typical relaxation time of the inter- 
acting gas, the existence of sound waves follows from the 
particle number conservation and Euler's equation (New- 
ton's equation for liquids) in a linear approximation^. In 
addition we have to assume here that these sound waves 
persist in the collision-less limit u ^> T~ elax , i.e. go over 
into acoustic phonons. For completeness I remark that 
besides of ordinary (first) sound in superfluids there is 
second sound which propagates at T — with the veloc- 
ity C2 = c/y/H. Second sound is a hydrodynamic mode 
which docs not exist in the collision- less limit 4 . 



II. T=0 

We begin with the properties of the ground state. In 
the non-interacting case all particle are in the state of 
lowest energy eo = and (^ ,T ^f') = uq = n where n = 
N/V. Here N and V denote the total number of particles 
and the volume, respectively. Pressure and the chemical 
potential vanish in the ground state. 

Phonons. As discussed already I will now assume that 
the excitations above the ground state can still be charac- 
terized by the momentum p (an assumption which par- 
allels Landaus picture of a Fermi fluid^) but that the 
dispersion relation e(p) is changed with respect to eo(p) 
such that 



qpi 

„2 



2m 



P^Po 



P > Po- 



(3) 
(4) 



The sound velocity c and the cross-over momentum po 
have to be determined self-consistently. Clearly at p = po 



cpo 



2m ' 



i.e. po ~ 2mc. 



(5) 
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I ignore here and in the following all factors of order unity 
since their determination is beyond the accuracy of our 
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considerations. It is also evident that for vanishing in- 
teraction the linear branch ([3]) of the dispersion relation 
should disappear. As long as the energy p 2 /(2m) of the 
free bosons is much larger than vqti (compare ([T])) the 
effect of the interaction is weak and the free boson dis- 
persion relation ([4]) will essentially remain unchanged. 
In the opposite case a linear dispersion relation is ex- 
pected. It is therefore obvious that a second relation for 
the determination of po is given byii 



2m 



v n 



mc 



(6) 



which gives with (JSJ) 



m m£;o 



Po ~ ^/v nm ~ — . (7) 
?0 



The relation ([7]) with C c (3) = 1 agrees for jj < ^ with 
the well-known resull££i£. Clearly the constant C c (d) (as 
well as other constants to come) remains undetermined 
in the present approach. Apparently £o ~ h/po plays 
the role of a fundamental length scale of the interacting 
system besides of the mean boson distance a = n~ x l d . £o 
is the smallest length scale on which collective behavior 
is still seen. In the superfmid phase it corresponds to the 
bare correlation or healing length. The assumption of 
weak interaction (which will be used everywhere below) 
can be written as 



a 



(d-2)/2 



< 1, 



-1/d 



(8) 



Zero point fluctuations. The degrees of freedom which 
are involved in the collective phonon excitations have a 
non- vanishing zero point energy since they correspond to 
an ensemble of harmonic oscillators. Hence they cannot 
contribute to the condensate. The depletion of the num- 
ber of particles No — Nno in the condensate is therefore 
given by 



riE>o = n - n a = 



J2 l = C n (d)n{~ 

p<Po 



(9) 



which agrees with the exact result in d = 3 dimensions 
if we choose Cjv(d = 3) = w 1.5 — . N denotes the 
number of particles in the condensate. Thus, contrary 
to the non-interacting case, the Heisenberg uncertainty 
relation in conjunction with the repulsive interaction en- 
forces that a finite fraction of the system is in excited 
states. 

The energy of the zero point fluctuations of the 
phonons follows analogously from: 



n e >oc po- 



rn 



P<PO 



The total ground state energy E consists then of the 
interaction of the particles in the condensate and the zero 



point fluctuations of the phonons 

Eo = -v N n + ]- ^ ^ c IpI 



P<PO 



vqNq 2 
2V 



I - Cruh ( - J 



(11) 



With C E (3) = yi^L s» 4.8 the last relation agrees 
with the exact result^. From (jTTj) follows the chemical 



potential 10 . 



M = 



dEo 
ON 



= v n 



1 + — 7-C E (d) - 1 

4 \ a 



The compressibility k — dn/ dpi is then 



1 + ^ CE (tY ,d - w 

8 V a J 



Finally we calculate the pressure 
dE Q 



P 



dV 



1 2 & 2 
2^° ~ „,„■ 



d+2 



(12) 



(13) 



(14) 



It is instructive to express the results ([9]) - (|14[) in terms 
of the healing length ©. 



III. T > 

Ideal gas. We begin again with a a brief review of the 
situation in the absence of any interaction 1 . Then at 
temperatures below the Bose condensation temperature 
of the ideal gas, 



T c ,o = Ct c - 



(15) 



whereCT c = 3, 31 11 , of the order N^a/Xxa) particles are 
outside of the condensate. Here At.2 = hj\JmT denotes 
the thermal de Broglie wave length of the non-relativistic 
Bose particles. The pressure follows as p ~ T/Xj. 2 etc., 
seei. 

Weak interaction, low T. For non-zero interaction 



Vono ~ mc 



^ rri (h 

T c ,o 



« T c .o (16) 



sets a second energy scale. As long as T <C mc 2 the rele- 
vant part of the excitation spectrum is essentially linear 
with a thermal de Broglie wave length 1 



he T cfi a? T Ci0 / b 



(d-2)/2 



(17) 



Since there is no phonon number conservation, Aj-j has 
the meaning of the average spacing between the phonons. 
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Since phonons behave as non-interacting relativistic par- 
ticle we obtain their entropy S by using the relation 1 
S = Cs(d)V/\^ ri for relativistic particles of energy T. 
The exact expression is C s (3) = 4tt 2 /90 w 0.44 i 2 -. 
Thus their energy E is given by 1 E — -^Cs{d)V 

and hence we get for the total free energy in the range 
T < mc 2 ~ T c . (b/a) d - 2 41 



T 



d+l 



(18) 

in agreement with the so-called Popov approximation^^. 
The finite temperature corrections are small, as has to 
be expected. Specific heat, internal energy, pressure etc. 
follow from that. 

Superfluid density. At zero temperature the whole 
Bose liquid is superfluid. At finite T however, part of 
it has normal viscous properties. To obtain the result 
for the normal and superfluid density, n n and n s , respec- 
tively, we follow the ideas of Landau and Khalatnikov^: 
We decompose the particle current density in a normal 
and a superfluid one 



n n Vr. 



(19) 



Here v s and v n denote the velocity of the normal and 
the superfluid component, respectively, in the laboratory 
frame. In the frame in which the superfluid is at rest the 
current is 



Jo 



,(v„ - v s ). 



(20) 



mj denotes the mass current density (i.e. the momen- 
tum density) of all particles which do not belong to the 
superfluid. This is the net momentum density resulting 
from all phonon degrees of freedom. It can be written as 

n n (v„ - v s ) = hr d I d d p —n p h {cp + p(v s - v„)) 
J m 

where n v h denote the phonon distribution function 
(which we don't have to know!). We expand the r.h.s. 
up to linear terms in 5~v = v s — v n , multiply the whole 
equation by <5v and perform the angular average. This 
gives 



1 

md 



= nCs{d) 



dep d mc 2 

, rfj-i , . d ( d - 2 ) 
T \ ( b\ — 5 — 



d d pcpn ph (p) 



S T 

I 

N mc 2 



(22) 



In the last steps we integrated by parts and expressed 
the resulting integral in terms of the energy and entropy, 
respectively. It is remarkable that the last relation is 
independent of dimension. Clearly the superfluid com- 
ponent has zero entropy. . 



Acknowledgment 

It is a pleasure to acknowledge helpful advice on the 
preparation of this paper by Jens O. Andersen, Sergey 
Artemenko, Andreas Engel, Axel Pelster, Aleksandra 
Petkovits, Zoran Ristivojevic, Friedmar Schiitze and 
Matthias Vojta. 



1 T. Nattermann, "A scaling approach to ideal quantum 
gases" Am. J. Phys 73, 349 (2005). 

2 In three dimensions (J5J) agrees with the textbook result 
in the limit of small wave vectors. For scattering in two 
dimensions see 

D.S. Petrov and G.V. Shlyapnikov, "Interatomic collisions 
in tightly confined Bos gas", Phys. Rev. A 64, 12706 
(2001). 

3 L.D. Landau and E.M. Lifshitz, "Fluid Dynamics", 2nd 
edition, Pergamon Press Oxford, §64 (1987). 

4 P.W. Anderson, "Basic Notions of Condensed Mat- 
ter Physics", Benjamin/Cummings Publishing Company, 
Menlo Park, California (1984). 

5 Most of the results for the weakly interacting Bose gas were 
originally obtained by T.D. Lee, K. Huang and C.N. Yang 
in "Eigenvalues and Eigenfunctions of a Bose System of 
Hard Spheres and Its Low- Temperature Properties" , Phys. 
Rev. 106, 1135-46 (1957). 

6 We refer here to the widespread textbook of L.D. Landau 
and E.M. Lifshitz, "Statistical Physics, Part 2", second 
edition, Pergamon Press Oxford (1981), eq. (25.11). 

7 More results on the weakly interacting Bose gas can be 



found in J.O. Andersen, "The Weakly Interacting Bose 
Gas", Rev. Mod. Phys. 76, 599-639 (2004). 

8 Referenced, eq. (25.19). 

9 Reference^, eq. (25.14). 

10 Reference^, eq. (25.15). 

11 L.D. Landau and E.M. Lifshitz, "Statistical Physics, Part 
1", third edition, Pergamon Press Oxford (1980), eq.(62.2). 

12 Reference^ 1 -, eq. (63.13). 

13 Referenced, eq. (27.10). 

14 Referenced, eq. (22.4). 

15 V.N. Popov, "Functional Integrals in Quantum Field The- 
ory and Statistical Physics", Reidel. Dortrecht (1983) 

16 see in particular I.M. Khalatnikov "An Introduction to the 
Theory of Superfluidity", W.A. Benjamin Inc, New York, 
Amsterdam (1965). 

17 Some of the formulas resemble relativistic counterparts , 
e.g. £ plays the role of the Compton wave length, mc 2 
that of the rest energy etc., but the crossover from non- 
relativistic to relativistic behavior is opposite to the con- 
ventional one: particles with low momenta behave as rela- 
tivistic, those with large momenta as non-relativistic ones. 



